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Abstract 

Shifts in species' traits across contrasting environments have the potential to influence ecosystem functioning. Plant 
communities on unusually harsh soils may have unique responses to environmental change, through the mediating role of 
functional plant traits. We conducted a field study comparing eight functional leaf traits of seventeen common species 
located on both serpentine and non-serpentine environments on Lesbos Island, In the eastern Mediterranean. We focused 
on species' adaptive strategies across the two contrasting environments and investigated the effect of trait variation on the 
robustness of core 'leaf economic' relationships across local environmental variability. Our results showed that the same 
species followed a conservative strategy on serpentine substrates and an exploitative strategy on non-serpentine ones, 
consistent with the leaf economic spectrum predictions. Although considerable species-specific trait variability emerged, 
the single-trait responses across contrasting environments were generally consistent. However, multivariate-tralt responses 
were diverse. Finally, we found that the strength of relationships between core 'leaf economic' traits altered across local 
environmental variability. Our results highlight the divergent trait evolution on serpentine and non-serpentine communities 
and reinforce other findings presenting species-specific responses to environmental variation. 
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Introduction 

Plants growing on special substrates (e.g. serpentine, limestone, 
gypsum, dolomite and shale) attract a lot of attention, not only due 
to their major contribution to global biodiversity but also because 
of their unique ecological character that may induce diverse 
community responses to environmental change [1]. Shifts in 
species traits across contrasting environments (e.g. productive vs 
rmproductive, polluted vs unpolluted) reveal the alternative 
strategies of plants for reproductive success and survival [2] and 
have the potential to influence ecosystem functioning (reviewed in 
[3,4]). For example, the harsh conditions of a substrate may limit 
the range of ecological strategies and thus filter the available 
species pools leading to communities dominated by species with 
similar functional traits [5,6] . Thus, dry environmental conditions 
may select for species with traits that allow them to use nutrient 
and water resources more conservatively [7]. However, although 
within each habitat abiotic environment leads to trait convergence 
by selecting similar trait values between coexisting species [8], 
niche differentiation leads to limiting similarity of trait values (trait 
divergence) [9]. In this context, plant communities on special 
substrates may be relatively responsive to changes (e.g. changes in 
rainfall, nitrogen deposition) due to their multiple limitations [10]. 
On the other hand, plants growing on special soils may be 



especially resistant to environmental changes due to their 
adaptations to harsh conditions [1 1]. Although shifts in functional 
traits across contrasting environments (e.g. wet-dry, productive- 
unproductive, etc.) may be expected [12] and have been captured 
by the major leaf economic dimensions [13,14,15], the variation of 
fundamental leaf traits across habitats and/ or within species is 
important to consider [16,17]. Moreover, although the relation- 
ships demonstrated by the leaf economic spectrum are robust at 
the global scale, the importance of trait variability on its core 
relationships is not well documented across locally contrasting 
environments [17]. 

Serpentine substrates are a well known example of a harsh 
environment for plants [18,19,20] and constitute efficient model 
systems for investigating variation on plant functional traits. 
Furthermore, serpentine ecosystems are important reservoirs for 
biodiversity as their flora includes a high number of rare and 
endemic species that present morphological and physiological 
adaptations to extreme conditions [18]. In a recent study, 
Californian serpentine grasslands showed greater resistance to 
environmental (climatic) fluctuation relative to non-serpentine 
[21], due to the presence of species with enhanced stress-tolerance 
traits (e.g. slow growth-rate, low height, low specific leaf area, high 
root/shoot biomass quotient; [22]). Mechanisms like abiotic stress 
and patchiness may also explain the greater temporal stability 



PLOS ONE I www.plosone.org 



1 



May 2014 | Volume 9 | Issue 5 | e96034 



Table 1. Leaf trait abbreviations and units. 
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(greater resistance to environmental fluctuations) of plant com- 
munities estabKshed in harsh environments [21,22,23,24,25]. 

Serpentine plant communities on Lesbos Island (eastern 
Mediterranean) have also shown higher short-term temporal 
stability in terms of species composition, relative to non-serpentine 
ones [26]. Although this may suggest a possible conservative 
response of these communities to environmental fluctuations (e.g. 
climate change [23,24]), it is not yet known if the higher short- 
term stability corresponds to species traits associated with efficient 
resource conservation. In this study we focus on species' adaptive 
strategies across serpentine and non-serpentine habitats on Lesbos 
and test across local contrasting environments the effect of trait 
variation on: a) the predictability of trait responses and b) the 
repeatability of relationships between core 'leaf economic' traits. 
Specifically, the following three questions are addressed: (1) Do 
species occurring on both serpentine and non-serpentine substrates 
present different adaptive strategies in response to different 
substrate types? If there is significant species leaf trait differenti- 
ation between the two contrasting substrates, will species occurring 
on serpentine substrates tend to have traits that allow them 
efficient resource conservation and species occurring on non- 
serpentine substrates tend to acquire resources rapidly? (2) Is there 
a repeatable ranking of species based on their leaf traits across 
different substrates? And finally, (3) if there is significant trait 
variation in response to substrate differentiation, are the relation- 
ships between the traits of the leaf economic spectrum conserved 
across contrasting substrates? 

Materials and Methods 

Ethics statement 

No specific permits and/or approvals were required for the 
described study sites on Lesbos (Greece). Given that aU the 
localities selected for our sampling are owned and managed by the 
government and are not private property or protected, no specific 
permits were required. In addition, our field study did not include 
any endangered or protected plant species. All data included in 
this study are freely available upon request. 

Study sites 

The study was conducted between May and June 2008 at four 
sites (Vatera, Ampeliko, Olympos and Loutra) located in the 
central and south-eastern part of Lesbos. All sites were dominated 
by herbaceous vegetation. A serpentine and an adjacent non- 
serpentine locality were chosen for comparisons in each of the four 
sites. The serpentine localities were selected based on how weU 
they represented the altitudinal and geographic range of serpen- 
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Figure 1 . Response patterns of SLA, LDMC, LT, LL, LW, LNC, LCC and LPC to different substrate types for seventeen herbaceous 
species (solid circles). Spearman ranl< correlation coefficients are given: ***, P<0.0001;**, P<0.01; *, P<0.05; not significant. The dotted line 
represents the 1:1 line. The bi-directional bars represent the standard error of means. Abbreviations are given in Table 1. 
doi:1 0.1 371 /journal.pone.0096034.g001 



tine habitats on tlie island and on the availability of accessible and 
adjacent non-serpentine areas within close proximity (0.6-7 km) 
and similar disturbance history and climatic conditions. AU 
selected non-serpentine localities were located on alluvial plains. 
A detailed description of the sites is available in Adamidis et al. 
[26] and in Kazakou et al. [27]. 

Leaf trait measurements 

We selected seventeen herbaceous species (Table SI) that were 
sufiiciently common and broadly distributed, all located in the two 
contrasting substrates. In order to capture the trait syndromes and 
thus the ecological strategies of these species, eight leaf traits 
(Table 1) were measured on the youngest fully expanded leaves 
with 10 replicates per species using standardized procedures [28]. 
Leaf length (LL) was measured as the distance from the leaf tip to 
the point of attachment with the stalk. Leaf width (LW) was 
measured as the diameter of the maximum imaginary circle fitted 
within the leaf [29]. Specific leaf area (SLA) was calculated as the 
ratio of the water-saturated leaf area to the leaf dry mass. Leaf dry 
matter content (LDMC) was determined as the ratio of leaf dry 
mass to water-saturated fresh mass. Leaf thickness (LT) was 
estimated by the (SLA x LDMC)"' product, or the water- 
saturated leaf fresh mass to leaf area ratio [30]. Leaf material was 
pooled to produce three batches that were ground separately and 
their leaf nitrogen (LNC) and leaf carbon (LCC) content were 
measured using an elemental analyser (Carlo Erba Instruments, 
model EA 1108, Milan, Italy). Data for leaf phosphorus 
concentration were obtained from Kazakou et al. [27]. 



Data analysis 

Where necessary, leaf traits were transformed to their natural 
logarithms before analysis to improve normality and homoscedas- 
ticity. Analyses of variance assuming species and substrate as fixed 
factors and locality as a random factor were used to explore the 
effects of species and substrate on all measured leaf traits. In the 
case of LPC, locality was not included in the ANOVA due to the 
sampling design and thus a two way ANOVA with both species 
and substrate as fixed factors was conducted. Spearman rank 
correlation coefiicients were calculated on untransformed leaf trait 
values to test whether the rank of species responses for each given 
trait was conserved across substrate types. For the evaluation of 
these bivariate correlations the species mean trait values were used 
for each substrate type (e.g. LDMC on serpentine substrates vs. 
non-serpentine). This analysis assesses the consistency of species 
ranking across different substrates based on species single-trait 
responses, i.e the species hierarchy in terms of mean trait values. 
To examine whether species shifts in response to different 
substrates were consistent across the multivariate-trait space a 
PC A ordination of the z-transformed species mean trait values was 
used for each substrate type. This analysis reduced our dimen- 
sionality in two axes (74% and 76% of the total variance explained 
for septenine and non-serpentine substrates respectively) and thus 
the position of each variable in two cUmensions. Each PCA 
produced a set of two-dimensional coordinates for each species-by- 
substrate combination, describing the position of each species on 
the multivariate-trait space for each substrate. In order to 
graphically represent the shifts of species leaf traits across 
multivariate trait space in response to substrate type differentia- 
tion, the coordinates of each species-by-substrate combination 




Axis 1 



Figure 2. Shifts of species leaf traits across multivariate trait space in response to substrate type differentiation. The grey dashed 
zero-centered vectors represent the species' shifts connecting the position of each species in leaf trait space across the two different substrate types. 
The black vector represents the mean shift of all species. 
doi:1 0.1 371/journal.pone.0096034.g002 
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Figure 3. Variation on relationsiiips between traits of tKie leaf economic spectrum across substrate type differentiation. Non- 
serpentine substrates are represented by solid circles and black lines, serpentine substrates are represented by open circles and grey lines, while the 
dashed lines represent the models that coincide. 
doi:1 0.1 371/journal.pone.0096034.g003 



were subtracted from each species coordinates on the serpentine 
substrate. In this way, the coordinates of each species on the 
multivariate-trait space for serpentine substrates are represented 
by the origin of the axes. Pearson correlation coefiicients were used 
for the investigation of among-trait relationships within each 
substrate type. 

Standardized major axis (SMA) analysis was used to examine 
whether relationships between core 'leaf economic' traits were 
conserved across different substrates. SLA-LNC, SLA-LPC, LPC- 
LNC and SLA-LDMC relationships were quantified and com- 
pared across different substrates. SMA was used to test whether 
these bivariate relationships differed in slope and when no 
significant differentiation emerged, differences in intercept (eleva- 
tion) or position along a common slope were tested. All the 
previous statistical analyses were carried out using the R statistical 
platform [31]. 

Results 

Effects of substrate type and species on leaf traits 

The ANOVA revealed significant effect of 'substrate' on all 
measured leaf traits (P<0.05; Table 2) except LNC and LPC {P> 
0.05). On average, SLA, LL and LW were respectively 20%, 



15.4% and 30"/o higher for individuals occurring on non- 
serpentine substrates while LDMC and LT were respectively 
7.8% and 22% higher for individuals occurring on serpentine 
substrates (Table 2). There were also significant 'substrate x 
species' (P<0.05; Table 2) interactions for all traits (apart from 
LCC), indicating that differences among species and localities 
affected the response of leaf traits to substrate types. 

Consistency of species' single- and multivariate- trait 
responses across different substrates 

Although the significant 'substrate x species' interactions denote 
that not all species respond similarly to substrate type transition, 
further analysis revealed significant across-species-correlations in 
mean values of leaf traits between the two contrasting substrates. 
Species presented a consisting ranking between serpentine and 
non-serpentine substrates for all leaf traits except for SLA, LNC 
and LPC (P>0.05) (Figure 1). 

On the other hand, the PC A of the z-transformed species mean 
trait values showed significant variation in species shifts across 
multivariate-trait space in response to substrate differentiation. 
The transition from serpentine to non-serpentine substrates caused 
widely diverse shifts in species' leaf trait values (Figure 2) that 



PLOS ONE I www.plosone.org 



6 



May 2014 I Volume 9 | Issue 5 | e96034 



Adaptive Strategies across Contrasted Environments 



resulted in a non-significant mean shift along the two dimensions 
(mean ±95% CI; Axis 1: 0.0059±1.03; Axis 2: 0.0059±0.82). 

Correlations among leaf traits and their consistency 
across different substrates 

SLA was negatively correlated with both LDMC and LT and 
positively correlated with LL, LW, LNC and LPC in the two 
substrate types (Table 3). SLA and LCC were positively correlated 
only on non-serpentine substrates (Table 3). LDMC was negatively 
correlated with LT and LPC and positively correlated with LL on 
both substrates. On serpentine substrates a negative correlation 
between LDMC and LW emerged. LL was positively correlated 
with LPC on serpentine substrates and negatively correlated with 
LT and LW on both substrates. LW was positively correlated with 
LNC on serpentine substrates and with LPC on non-serpentine 
substrates. Finally LNC was positively correlated with LCC on 
non-serpentine substrates and with LPC on both substrates 
(Table 3). 

Standardized major axis (SMA) analysis was used to examine 
the consistency of four relationships between core 'leaf economic' 
traits across different substrates. Significant differentiations in 
slopes across different substrates emerged for the relationships SLA 
vs LNC, LPC vs LNC and SLA vs LDMC (Table 4; Figure 3). The 
three relationships were significant on both non-serpentine and 
serpentine substrates (Table 4). The models describing the 
relationship SLA vs LPC coincided across different substrates 
(Table 4; Figure 3). The relationship SLA vs LPC was significant 
across serpentine and non-serpentine substrates (Table 4). 

Discussion 

All measured leaf traits, except LNC, LCC and LPC, varied 
E significandy across the different substrates. SLA, LL and LW 

presented on average higher values on non-serpentine substrates 
while species from serpentine substrates showed higher values of 
LDMC and LT. Species trait values from non-serpentine 
substrates are associated with a more exploitative strategy. 
Exploitative species tend to acquire resources rapidly by presenting 
high values of SLA (low values of LDMC) along with high relative 
growth and photosynthesis rates [15]. In addition, it is well known 
that leaf sizes (leaf length and leaf width) are smaller in drier and 
nutrient poor environments [7,32] such as serpentine habitats. On 
the other hand, species from serpentine substrates presented a 
conservative strategy, investing more resources to structural 
compounds and thus presenting higher values of LDMC (low 
values of SLA) and leaf thickness (denser leaves). Considering that 
leaf thickness has been associated with water storage processes 
[33], the higher values of this trait on species inhabiting serpentine 
substrates may relate to the low water-holding capacity of these 
substrates [34]. Hence with respect to our first question, the 
transition between the two contrasting environments is associated 
with changes in traits which are consistent with the leaf economic 
trade-off [14,15,35] and highlight differences in functional 
strategies followed by the same species at different environments. 

In our study, values of LCC did not significantly differ between 
serpentine and non-serpentine substrates despite the fact that high 
values of this trait represent investments in structural strength [36]. 
Grassein et al. [37] suggest that species' strategies are associated 
not only with plant traits but also with trait plasticity, demon- 
= g ^ strating that conservative species exhibit constant values of a 
"I ^ structural trait across a resource availability gradient in contrast to 
I c 5 a significant range of values presented for this trait by the 
< o ;5 exploitative species. In our case, it is possible that higher trait 
plasticity of species from non-serpentine substrates in LCC is 
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responsible for this lack of differentiation across contrasting 
substrates. In addition, both LNC and LPC were unaffected by 
the substrate type. It appears that the effect of any potential 
differentiation in nutrient availability between serpentine and non- 
serpentine substrates is diluted due to the overall low nutrient 
availability that characterizes Mediterranean ecosystems 
[27,36,38,39]. Navas et al. [36] also found no significant variation 
in leaf nitrogen and phosphorus concentrations across Mediterra- 
nean successional stages and attributed this lack of differentiation 
to the low soil nutrient availability. 

The significant species by substrate interactions that emerged 
for all leaf traits (except LCC; Table 2) indicate the variation in 
species responses across the different environments. However, the 
strength of this interaction's effect depends on the trait under 
consideration. These results are consistent with other studies 
demonstrating differentiation in species responses to environmen- 
tal variation [40] and may reinforce studies suggesting that within- 
species trait variabihty is not only considerable but also species- 
and trait-specific [17,41]. 

Despite the fact that species differed in their trait responses to 
substrate types, it seems that at the single-trait level, species 
responded in the same direction. All the single-trait species' 
correlations between the two contrasting environments, except for 
SLA, LNC and LPC, were significant (Figure 1). In other words, 
with respect to our second question, the species' single-trait 
responses were generally consistent, conserving the species ranking 
across the contrasting environments. Consistent rankings of 
species' trait responses have also been demonstrated across several 
spatial, temporal, environmental and climatic gradients 
[42,43,44,4,5]. On the other hand, the transition from serpentine 
to non-serpentine substrates caused widely diverse species' 
multivariate-trait responses that tended to cancel each other out 
and result in a non-significant mean species shift (Figure 2). Thus, 
the transition between the two contrasting environments forced 
the species' traits to respond idiosyncratically. This result 
reinforces the generalization of the pattern found by Wright and 
Sutton-Grier [40], who also demonstratixi uncoordinated trait 
variation while studying a different set of species under control 
conditions. The coordinated leaf trait variation described among 
global vegetation [15,46,47,48,49] is not supported in our local 
scale study between serpentine and non-serpentine habitats. 

The significant bivariate correlations that emergc-d hc-tween 
several leaf traits within both serpentine and non-serpentine 
substrates followed similar patterns across the two contrasting 
environments (Table 3). The relationships between traits of the leaf 
economic spectrum were significant on both substrate types 
(Table 4). However, the SMA analysis revealed significant slope 
differentiation for the relationships SLA vs LNC, LPC vs LNC and 
SLA vs LDMC across the different substrates. Our results indicate 
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